Pathways to double ionization of atoms in strong fields 
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We discuss the final stages of double ionization of atoms in a strong linearly polarized laser field 
within a classical model. We propose that all trajectories leading to non-sequential double ionization 
pass close to a saddle in phase space which we identify and characterize. The saddle lies in a two 
degree of freedom subspace of symmetrically escaping electrons. The distribution of longitudinal 
momenta of ions as calculated within the subspace shows the double hump structure observed in 
experiments. Including a symmetric bending mode of the electrons allows us to reproduce the 
transverse ion momenta. We discuss also a path to sequential ionization and show that it does not 
lead to the observed momentum distributions. 
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I. INTRODUCTION 

Present day lasers are powerful enough to ionize several 
electrons from an atom. The electrons can be removed 
one by one in a sequential process or all at once in a non- 
sequential process. Independent electron models give ion- 
ization rates that are much smaller than experimentally 
observed indicating that interactions between elec- 
trons are important. A series of most recent experiments 
has added the observation that also the final state of 
the electrons is dominated by the interactions: the to- 
tal momentum of the electrons is aligned along the field 
axis and the joint distribution of parallel momenta 
for the two electrons, in the double ionization experi- 
ment has pronounced maxima along the diagonal, 
showing that the electrons typically come with the same 
momenta. These observations have been reproduced in 
numerical simulations with varying approximations and 
simplifying assumptions But given the complexity 

of the analysis that is required the essential elements are 
difficult to identify. As a step towards a better under- 
standing we discuss here the pathways to double ioniza- 
tion within a classical model for electrons in a combined 
Coulomb and external field. 

Our aim is not to describe the full ionization process 
all the way from the ground state to the final, multiply 
ionized state. According to the currently accepted mod- 
els f§|~|l3l,|6|,f7f , the whole process of multiphoton multiple 
ionization can naturally be divided into two steps: in the 
first step a compound state of highly excited electrons 
close to the nucleus is formed and in the second step sev- 
eral electrons can escape from this compound state to 
produce the multiply ionized final state. We focus on the 
last step, the escape of two or more electrons from the 
highly excited compound state close to the nucleus. 

The formation of the intermediate compound state is 
suggested by the rescattering model J||J for strong field 
multiple ionization. According to this model the en- 
hanced cross section for multiple ionization is due to a 
rescattering of one electron that is temporarily ionized 



and accelerated by the field before it returns to the nu- 
cleus when the field reverses. During the collision energy 
is transfered to other electrons, but all of this happens 
close to the core, where the dynamics of the electrons is 
fast and the interactions are strong and non-integrable. 
As a consequence, details of the initial preparation pro- 
cess are lost. Moreover, the decay of this state is also 
quick. The compound state is thus a short lived, highly 
unstable complex that separates the first half of the ex- 
citation process, whose main contribution is the build up 
of energy in the complex, from the second half, where the 
decay mode is determined. 

The compound state has several decay paths: single 
ionization when only one electron escapes to infinity, dou- 
ble or multiple ionization with two or more electrons es- 
caping to infinity, and the case of a single electron that es- 
capes from the neighborhood of the nucleus but is rescat- 
tered to the next field reversal: in the latter case the 
whole process repeats itself, another compound state is 
formed and the decay path has to be selected anew. 

We discuss here the further evolution after the forma- 
tion of the compound state. To be specific, we will focus 
on the escape of two electrons in the following, but the 
arguments can easily be extended to the removal of more 
than two electrons. The main aim is to identify the chan- 
nels that lead to double escape and to study their signa- 
tures in the distribution of electron and ion momenta. 
Our analysis is purely classical. Given the highly excited 
complex from which we start and the multiphoton nature 
of the process this seems a reasonable point of entry to 
the final stage of the ionization process. Our analysis is 
very similar to Wanniers approach |L1|-|T7|] to double ion- 
ization by electron impact. There main difference is that 
in the present case one has to take into account also the 
external field. Brief summaries of some aspects of this 
model have been presented in Jl^,|l{| . 

The pathways to double ionization are discussed in sec- 
tion 2. The dynamics in the Ci v and C v subspaces in- 
cluding the effective potential and sample trajectories is 
discussed in section 3. A key element of our argument is 
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the identification of a saddle in phase space near which 
trajectories leading to non-sequential double ionization 
have to pass; its properties and stability are discussed in 
section 4. The distributions of electron momenta within 
the Civ and C v subspaces are analyzed in section 5. The 
dynamics outside these symmetry spaces and sequential 
ionization processes are discussed in section 6. We con- 
clude with a summary of the model in section 7. 



II. PATHWAYS TO DOUBLE IONIZATION 

As described in the introduction we assume an ini- 
tial state of two highly excited electrons near the nucleus 
which then decays to cither single ionization or double 
ionization. During this process the linearly polarized 
laser field is always on. Therefore, the Hamiltonian con- 
sists of three parts, 



H = T + V l + Vi 
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the kinetic energy of the electrons, 



9 2 

T= Pi K 
2 2 



(1) 



(2) 



the potential energies associated with the interaction 
with the nucleus and the field (polarized along z-axis), 

v * = -]t\~W\ +Fit)zi + mz2 (3) 

and the repulsion between the electrons, 
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Vl2 



r 2 - ri 



(4) 



The electric field strength F(t) has an oscillatory com- 
ponent times the envelope from the pulse; the discussion 
applies to general F(t), and specific choices have to be 
made for the numerical simulations only. 

Once the electrons leave the atom the repulsion pushes 
them apart and becomes weaker the larger the separa- 
tion. Thus in the asymptotic state after ionization and 
after the pulse is turned off, repulsion is minimized. In 
order to identify the effects of the electron repulsion on 
the full process it is instructive to consider the double 
ionization events without repulsion first. This will be 
presented in the next section. The pathways with elec- 
tron repulsion included will be presented thereafter. 



A. Without electron repulsion 

The Hamiltonian for two electrons without electron re- 
pulsion splits into two independent Hamiltonians for each 
electron. In view of the fast motion of the electrons close 
to the nucleus we will frequently use an adiabatic as- 
sumption and discuss motion of the electrons in a field 



with fixed amplitude F. Note, however, that all simula- 
tions use the full time-dependent field and do not make 
use of this adiabatic assumption. With this assumption 
each electron moves in a constant electric field, one of the 
few non-trivial integrable problems. The initial energy 
and the other constants of motion of each electron are 
fixed by the initial conditions (in the compound state) 
and do not change. Double ionization can thus occur 
only if both electrons individually have enough energy 
to ionize (and have the other constants of motion so as 
to allow ionization). The threshold for ionization is set 
by the field strength: if the field is non-zero a Stark 
saddle forms and the total energy has to be above the 
Stark saddle. F or a co nstant field strength F the saddle 
lies at \zp\ = y/2/\F\ and has a potential energy (single 
electron) Vp = —2y/2\F\. Therefore, double ionization 
is excluded if the total energy for both electrons is less 
than 2Vp. For E = 2Vp the only path leading to double 
ionization has both electrons with the same energy. For 
E > 2Vp double ionization becomes possible even with 
slight asymmetries in the energy distribution. 

In the full system of two electrons without repulsion 
but in a time-dependent field integrability is lost but sep- 
arability is still preserved. As discussed in the introduc- 
tion the initial state is a compound state with negative 
energy for each electron. If that energy is above the Stark 
saddle and if the motion is directed towards it the elec- 
trons can cross and run away from the nucleus. Once 
they cross the saddle they have to gain energy to escape 
from the Coulomb attraction when the field is turned off. 
This happens essentially by running down the potential 
energy slope on the other side of the saddle while the field 
is still on. This mechanism of energy gain is the same as 
in the interacting electron case. 



B. With electron repulsion 

With the electron repulsion included the common 
Stark saddle at zp is no longer accessible since the elec- 
trons cannot sit on top of each other. The best that can 
be achieved is a symmetric arrangement of both elec- 
trons in the same distance from the nucleus and sym- 
metric with respect to the field axis. If the distances 
are not the same, electron repulsion will push the elec- 
tron that is further out away from the nucleus and thus 
help towards ionization, but the one further in has to 
face not only Coulomb attraction but also the repulsion 
from the one further out. Thus repulsion acts so as to 
amplify differences in energy in this configuration. The 
configuration that is singled out is a symmetric one, with 
both electrons moving at the same distance from the nu- 
cleus. Deviations from this configuration will be ampli- 
fied sufficiently so that non-sequential double ionization 
is suppressed and only single ionization takes place. A 
remaining electron can, however, be still ionized if its en- 
ergy is higher than the saddle energy for a single electron 
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atom and, in the adiabatic picture, the other constants of 
motions allow to do so. In such sequential ionization cor- 
relations between escaping electrons are strongly weak- 
ened, that is the final momenta of the electrons along the 
polarization axis can be either parallel or anti-parallel. 



III. DYNAMICS IN THE C 2 v AND CV 
SUBSPACES 



the pulse envelope of duration T^. 

Setting r = 0, p r = 0, z = and p z = we define the 
C v symmetry subspace. The Hamiltonian is then reduced 
to 



where 



H = 



P% 



P% 



(9) 



A two electron atom illuminated by a linearly polar- 
ized electromagnetic wave possesses some symmetry sub- 
spaces. The simplest Civ symmetric configuration corre- 
sponds to both electrons moving in a plane which con- 
tains the field polarization axis and with positions and 
momenta symmetric with respect to this axis. The elec- 
trons put in such a configuration never leave it because 
there is no force which can pull them out of the subspace. 
The symmetry subspace can be enlarged. That is, with 
additional bending motion of the electrons with respect 
to the field axis the symmetry subspace is C v . 

The symmetric configurations become more apparent 
in suitably chosen coordinates. We apply the canonical 
transformation 



R = (Pl +P2)/2, PB.=P Pl +Pp 2 
r =(Pl-P2)/2, Pr=P Pl -Pp 2 

Z=(zi + z 2 )/2, pz=p Zl +Pz 2 
z= (zi- z 2 )/2, Pz=p zl ~Pz 2 , 
ip = epi + <p 2 , L = (p tfil +p ip2 )/2 
<p = tp 1 - ip 2 ,P4, = (jVi ~Pvi)l 2 



V = — 



(5) 



where (pi,Zi,(pt) are cylindrical coordinates of the elec- 
trons, labeled i = 1 and 2. For double ionization in lin- 
early polarized laser field the total angular momentum 
projection on the polarization axis is conserved. The ex- 
periments begin with atoms in the ground state, thus, for 
the field directed along the z-axis we may choose L = 0. 
Then the Hamiltonian of the system reads 



H 



Pr+pI+Pz+pI 



P\ 



2{R-r) 
with the potential energy 



2{R + r) 2 
2 +V(R,r,Z,z,<p,t) 



(0) 



V = -- 



V(i? + r) 2 + (Z + z) 2 y/(R-r) 2 + (Z - zf 
1 



^2R 2 - 2{R 2 
-2ZF(t), 



r 2 ) cos (j> + 2r 2 + 4z 2 



VR 2 + Z 2 i?^2(l-cos< 



+ 2ZF(t). (10) 



The potential (^) is symmetric in r (for z = 0) and in z 
(for r = 0), so that the derivatives with respect to r and z 
vanish: once the electrons are in the symmetry subspace 
r = z = and p r = p z — 0, they cannot leave it. 





where the field is given by F(t) = Ff(i) cos(ujt + 9) with 
F, lu and 9 the peak amplitude, frequency and initial 
phase of the field, respectively, and with 



f(t) =sm 2 (Trt/T d ) 



(7) FIG. 1. Sections through equipotential surfaces of the adi- 
abatic potential Eq. (^) for fixed time t and for <j> = n (a) 
and (j> = 7r/4 (b). Panel (a) corresponds also to the poten- 
tial Eq. ( |12[ ) in the Ci v symmetric subspace; the saddle moves 
along the dashed line when the electric field points in the pos- 
itive Z-direction and along a second obtained by reflection on 
(8) Z — during the other half of the field cycle. 
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The further constraint <fi — tt and p$ = leads to the 
Civ symmetry subspace 



h= Pr + P% +v{r ^ 



with potential 



V 



4 1 
Vi? 2 + Z 2 + 2R 



2ZF{t). 



(11) 



(12) 



Let us begin with an analysis of the motion in the Ci v 
subspace. The electrons move in a plane and their posi- 
tions {pi — R, zi = Z) in that plane and their momenta 
(j?Pi = P-r/ 2 , p Zi = Pz/2) are the same. 

The adiabatic potential (Q) for fixed time correspond- 
ing to the maximal field amplitude F — 0.137 a.u., i.e. 
an intensity of 6.6 • 10 14 W/cm 2 , is shown in Fig. [l]a. 
The saddle is located along the line Zg = r$ cos 0s and 
Rs = r s sin 9s with 9 s = 7r/6 or 5ir/6 and at a distance 



r 2 =V3/\F(t)\. 



The energy of the saddle is 



V s = -6J\F(t)\/V3. 



(13) 



(14) 



For the above mentioned field the saddle has an energy 
of V s = -1.69 a.u.. 

If we switch off the repulsion between the electrons the 
saddle will move onto the Z-axis, i.e. both electrons are 
allowed to escape symmetrically on top of each other. 
With repulsion the saddle splits into two symmetrically 
related ones and moves away from the Z-axis. 

A typical trajectory within the symmetric configura- 
tion for uj — 0.057 a.u. (800 nm) is shown in Fig. |^. 
During the ramping of the field the electronic motion is 
little influenced by the electric field, but during the third 
half cycle of the field the saddle is close enough to the 
electron orbits and ionization takes place. Once on the 
other side of the saddle, the electrons rapidly gain en- 
ergy. The saddle thus provides a kind of transition state 
p0| , pl| for the correlated double ionization process: once 
the electrons cross it, they are accelerated by the field 
and pulled further away, making a return rather unlikely. 
Moreover, the electrons can acquire the missing energy 
so that both can escape to infinity even when the field 
vanishes. Note that before double ionization occurs the 
effect of the field on the electrons is minimal, supporting 
the adiabatic assumption. 

In the experiments (2]-[j| ion momenta both parallel 
and transverse to the field are measured. For uj = 
0.057 a.u. momentum transfer by photons is negligible, 
so the ion momentum reflects the sum of the momenta 
of the emitted electrons, Pi + P2 = ~Vion- Symmet- 
ric motion in the Ci v subspace takes place in a plane 
and consequently the total transverse momentum of the 
electrons is zero. In the C v subspace the electrons are 



permitted to leave the plane, and this is a minimal ex- 
tension necessary to give obtain non- vanishing transverse 
momenta. 
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FIG. 2. A typical trajectory in the Ci v symmetric sub- 
space with E — —1.3 a.u., F = 0.137 a.u. and the pulse 
duration — 4 x 2ir/uj. Panel (a) distance of the electrons 
to the nucleus. The dashed line indicates the distance of the 
saddle. Note that before the double ionization occurs the ef- 
fect of the field on the electrons is minimal, supporting the 
adiabatic assumption. Panel (b) energy of the electrons. Note 
that the initial state has negative total energy and cannot lead 
to double ionization. The energy increases once the electrons 
have escaped from the nucleus far enough so that acceleration 
by the electric field dominates. 

The sections through equipotential surfaces of Eq. ( |To| ) 
for cj) — 7T and <f> — it /A are shown in Fig. |l|. Increasing or 
decreasing 4> from tt results in greater repulsion energy 
between the electrons. So the electrons can escape, with 
the coordinate </> ^ 7r, only if the energy is greater than 
the saddle energy, Eq. (14). 

We discuss here double ionization of atoms but the 
analysis of a symmetric escape can be easily extended 
to a multiple ionization process. Especially for TV elec- 
trons the symmetric subspace corresponding to N parti- 
cles symmetrically distributed in the plane perpendicular 
to the polarization axis is Cn v J19|. 



IV. THE SADDLE 

Within the adiabatic assumption and for a fixed exter- 
nal field strength non-sequential double ionization hinges 
on the crossing of the saddle in the Ci v subspace. For 
the dynamics within the subspace this is obvious from the 
potential. For the motion in the full six-dimensional con- 
figuration space this is less clear. One possibility to gain 
insight locally near the saddle is a harmonic analysis and 
a determination of the frequencies of small deviations. 
Within the Ci v symmetric subspace there is one hyper- 
bolic mode corresponding to motion across the saddle 
(th e 'reaction coordinate', in chemical physics parlance 
J,|l|) and a stable one corresponding to perpendicular 
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variations. In the full space we expect at least one ad- 
ditional unstable one, corresponding to the amplification 
of energy differences mentioned before in section (IIB). 
The analysis in this section is for fixed field strength, for 
electron dynamics in a constant external field, justified 
by the adiabatic reasoning. 

The second order variations of the potential (Q) around 
the saddle point results in 



H 



P%+Pz 



4 
3^3 



zr s 



t-(r - r s ){z - z s ) - ^r(z - z s y 



+ 



r s 

pI + p 2 z 



2r| 



4Rc 



3V3 



tt) 2 + V s . 



2r|' 



(15) 



The (/)-degree of freedom corresponds to a bending mo- 
tion of the electrons against each other and is stable on 
account of the repulsive nature of the Coulomb force. 
This degree of freedom is also the one that comes in by 
going from the C-^v subspace to the C v subspace. 

Diagonalization in the (R — Rs, Z — Zs) space reveals 
one stable and one unstable mode. The latter corre- 
sponds to the reaction coordinate and its Lyapunov expo- 
nent is /i w I.2IF 3 / 4 . Similar analysis in the (r, z) space 
yields another stable and unstable mode with Lyapunov 
exponent v w 1.57F 3 / 4 . The direction of the unstable 
mode is (w r ,w z ) ~ (0.43u>,u>) and it corresponds to the 
situation when one electron escapes and the other one is 
turned back to the nucleus. That is, with positive and in- 
creasing w the first electron moves away from the nucleus, 
i.e. pi = Rs + w r and Z\ = Zs + w z grow [see (§], while 
the other one returns to the nucleus, i.e. p\ = Rs — w r 
and z\ = Zs — w z decrease. 

All in all there are three stable modes and two un- 
stable ones. Any energy contained in the stable modes 
is preserved and cannot be transferred to kinetic energy 
along the reaction coordinates. Whether single or dou- 
ble ionization occurs is thus determined solely by the 
energy and initial conditions in the two hyperbolic sub- 
spaces. For energy equal to the saddle energy only a 
trajectory within the Civ symmetry subspace leads to 
non-sequential escape - any deviation from the subspace 
are growing faster than the escape along the reaction 
coordinate since v > p. For energy higher than the sad- 
dle some deviations from the symmetry subspace are al- 
lowed. In particular, following Wanniers lead and Rosts 
generalization pj-p"7| it is possible to estimate the crit- 
ical behavior for the double ionization cross section at 
threshold. It is algebraic with the exponent given by the 
ratio of the positive Lyapunov exponents of the unstable 
modes. A detailed discussion of this is outside the main 
line of our arguments here and will be given elsewhere. 

It is instructive to actually calculate numbers for the 
Lyapunov exponents in the two directions. For lasers 



with the maximal field strength of F — 0.137 a.u. we 
find l/[i = 3.7 a.u. and 1/u = 2.8 a.u.. Compared to the 
period of the laser, 2tt/w = 110.2 a.u. this is rather fast, 
indicating that the crossing of the saddle and the separa- 
tion away from the double ionization manifold take place 
rather quickly. This justifies also our adiabatic analysis 
in this section. 



V. FINAL STATE MOMENTA DISTRIBUTION 
WITHIN THE SYMMETRIC SUBSPACES 



So far we have discussed the phase space features in an 
adiabatic approximation for fixed field strengths. Now 
we will use this to draw conclusions about the experi- 
mentally observed signatures, specifically about the dis- 
tributions of ion momenta in the final state. They can 
be calculated within the C v subspace by averaging over 
all initial conditions and all phases of the field. That this 
is possible is connected with the instability of the saddle: 
all trajectories leading to the non-sequential ionization 
have to pass sufficiently close to the saddle and the sym- 
metric subspace. It therefore is possible to estimate the 
behaviour near the subspace from the dynamics within 
the subspace. 



A. Parameters of the model 



The rescattering of an electron leads to a highly excited 
complex of total energy E which is one of the parame- 
ters of our model. The maximal energy a rescattering 
electron can bring in has been estimated to be about 



Z.llUp |§|9|], where U p is the ponderomotive energy of an 
electron. For the weakest field used in the experiment 
on double ionization of He atoms 0] this maximal en- 
ergy barely corresponds to the ionization energy of the 
He 1+ ions. We therefore assume in the following that 
the highly excited complex has a negative initial energy, 
E < 0. 

The absence of detailed knowledge of the structure of 
the initial compound nucleus suggests to average over the 
initial configurations. However, even for a negative en- 
ergy and fixed time it is difficult to define a microcanon- 
ical distribution of initial conditions for the Hamiltonian 
(^) since, for non-zero external field, the system is open. 
Therefore, we choose for the calculations initial condi- 
tions from the energy shell that also lie in the hypersur- 
face Z = 0. The results are not sensitive to a particular 
choice of the hypersurface but the one for Z = has the 
advantage that the dipole moment along the filed is zero 
and the choice of the initial conditions does not depend 
on the initial field phase. 

The second parameter, in addition to the energy, is the 
time to during the pulse Eq. (||) when the highly excited 
complex is formed. The rescattering event is not possible 
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at the beginning of the pulse, so one has to start simula- 
tions somewhere in the middle of the pulse. In Fig. |^ and 
H final distributions of ion momenta for the initial energy 
E = -0.1 a.u., field strength F = 0.137 a.u. and dif- 
ferent initial time to are shown. The distributions of the 
transverse momenta are almost the same but the ones for 
the parallel momenta differ. The latter reveal a double 
hump structure with widths sensitive to the initial time. 
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FIG. 3. Final distributions of ion momenta parallel, p par 
and transverse, ptrans to the field polarization axis for the ini- 
tial energy E = —0.1 a.u., peak field amplitude F = 0.137 a.u. 
and pulse duration Td = 20 x 2ir/u). Panels (a)-(b) cor- 
respond to the initial time to = 0.257d in the pulse dura- 
tion with the envelop f(t) — s'm 2 (irt/Td) while panels (c)-(d) 
to to = 0.75Td. Dashed lines are related to the estimates 
±2Ff(to)/u> = ±2.4 a.u.. Note that the distributions are es- 
sentially the same independently if one chooses to before or 
after the peak field value provided /(to) is the same. The 
results are based on integrations of about 8 ■ 10 4 trajectories. 
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FIG. 4. The same as in Fig. | but for to = 0AT d 
[panels (a)-(b)] and to = 0.6T d [panels (c)-(d)]. The width 
of the parallel momentum distribution can be estimated as 
±2Ff(t )/uj = ±4.3 a.u.. 



The maximum energy that can be acquired by a free 
electron in the field is 2U P . So, for parallel emission of 
two electrons, the maximal parallel ion momentum can 
be estimated as 



P P ar^2^W p = 2Ff(t)/uJ. 



(16) 



If we substitute in Eq. ( |l6| ) t = to we find values which 
correspond very well to the widths of the distributions in 
Fig. || and ||. The figures indicate also that the widths 
are the same independently if to is chosen before or after 
the peak of the pulse provided /(to) is the same. This 
implies that the dominant ionizations take place in the 
first field cycle after the complex is formed. Fitting the 
width of the calculated distribution to the experimental 
results allows one to estimate the moment in the pulse 
when majority of doubly ionized ions are created. 
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but for fixed t = 0.33T d 



-3 3 6 

Ppar (a-U.) 
FIG. 5. The same as in Fig. 
and different initial energy: E = —0.05 a.u. [panels (a)-(b)], 
E = —0.5 a.u. [panels (c)-(d)] and E = —1.5 a.u. [panels 
(e)-(f)]- 

Now we fix to and change the initial energy E; the 
results are shown in Fig. B. For the lowest energy, 
E = —1.5 a.u., the transverse momentum distribution 
is narrower than for higher energies. E = —1.5 a.u. is 
actually very close to the minimal saddle energy for the 
applied field, Vs — —1.69 a.u.. Thus the effect is natural 
as, close to the saddle energy, only trajectories near the 
Civ subspace can cross the saddle and those with <j) ^ it 
bounce back from the potential barrier and do not ionize. 

While the width of the distributions of the parallel 
momenta do not change significantly their shapes do, es- 
pecially for the initial energy close to the saddle. The 
electrons that cross the saddle when the energy is near 
Vs are slow and the combined interactions of the external 
and Coulomb fields shapes the distributions. This is dif- 
ferent from the high energy case: then shortly after the 
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electrons cross the saddle the interaction with the elec- 
tric field is stronger than the attraction to the nucleus 
and the distribution is mostly shaped by the laser field. 

The initial energy of the complex is the higher the 
higher intensity of the laser is and the larger the energy a 
rescattered electron can bring in ^|J9). From the depen- 
dence of the distributions on the initial energy we may 
conclude that in the experiment the shape of the distri- 
bution of the parallel ion momenta should change char- 
acter when the laser intensity increases. For the intensity 
at the threshold for non-sequential double ionization the 
distribution with single maximum around zero momen- 
tum is expected; for higher intensities the double hump 
structure should turn up. 

All numerical results have been obtained for initial con- 
ditions taken from the C v symmetry subspace. Our re- 
sults for the C% v subspace for parallel momentum dis- 
tributions are essentially the same. The transverse mo- 
menta of ions for the C2 V subspace are, however, zero 
because of the symmetry assumption. 

After this discussion of the two parameters (initial en- 
ergy and starting time of the integration) we can turn to 
comparisons with experimental observations. 



B. Comparison with experimental results 

Weber et al. Q carried out double ionization exper- 
iments with He atoms and measured the distributions 
of ion momenta. They applied infrared (800 nm) laser 
pulses of the duration 220 fs (measured on the half 
peak value) and with the peak intensities in the range 
(2.9 - 6.6) x 10 14 W/cm 2 . In Fig. § and g we show the 
experimental distributions and compare with those cal- 
culated in the C v subspace. The agreement is very good 
except for the parallel momentum distribution in Fig. 
where the calculated distribution possesses a much more 
pronounced minimum than in the experiment and the 
positions of the peaks do not exactly coincide with the 
experimental values. 

There are a few possible sources for these discrepan- 
cies. First, the pulse duration in the experiment was 
quite long, i.e. about 80 field cycles. The slow ramp- 
ing of the field in the experiment implies that the initial 
time to of the ionization is less well defined, i.e. there are 
contributions from some range of to- There are also con- 
tributions from different initial energies. Secondly, real 
ionizing trajectories do not live exactly in the symmetry 
subspace but close to it, leading to asymmetries and ad- 
ditional differences in the final momenta. And there are 
also possible contributions from sequential double ioniza- 
tion events (see below). 

Moshammer et al. |4[] performed experiments with Ne 
atoms for much shorter pulses, i.e. 30 fs and for radiation 
with similar wave length (795 nm) as the previous group. 
The comparison of our calculations with this experiment 
are presented in Fig. K Interactions of the two excited 



electrons with the other electrons are neglected in our 
model and the energy values used in the calculations are 
measured with respect to the threshold for the two elec- 
tron continuum (i.e. about 2.3 a.u.). The agreement is 
even better than for the case of He atoms. This is pre- 
sumably due to the much shorter pulse duration and the 
fasted ramping of the field, so that the time to when the 
majority of the ionization events happen is much better 
defined. 




FIG. 6. Panels (a)-(b): final ion momentum distributions 
measured in the experiment of double ionization of He atoms 
in the focus of 800 nm, 220 fs (i.e. about 80 x 2-k/ui) laser 
pulses at peak intensity of 2.9 x 10 14 W/cm 2 (i.e. for the field 
strength F = 0.091 a.u.) from |). Panels (c)-(d): the cor- 
responding distributions calculated in the C„ symmetry sub- 
space for the initial energy E = —0.6 a.u. and to = 0.33T<i 
where T d /2 = 80 x 2tt/lu, see Eq. §. 
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FIG. 7. Panels (a)-(b): the same as in the corresponding 
panels in Fig. ^ but for the peak intensity of 6.6 x 10 14 W/cm 2 
(i.e. for the field strength F = 0.137 a.u.). Panels (c)-(d): the 
same as in the corresponding panels in Fig. ^ but for the initial 
energy E — —0.4 a.u.. 
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FIG. 8. Panels (a)-(b): final ion momentum distributions 
measured in the experiment of double ionization of Ne atoms 
in the focus of 795 nm, 30 fs (i.e. about 11 x 2-k/ui) laser 
pulses at peak intensity of 13 x 10 14 W/cm 2 (i.e. for the 
field strength F = 0.192 a.u.) from [§. Panels (c)-(d): the 
corresponding distributions calculated in the C v symmetry 
subspace for the initial energy E = —0.3 a.u. and to = OATd 
where T d /2 = 11 x 2n/u, see Eq. |. 



VI. SEQUENTIAL DOUBLE IONIZATION 



Already from the experiment it is clear that double 
ionization is a rare process, e.g. outnumbered by single 
ionization by about 10 4 : 1 for He atoms and field inten- 
sity 2.9 • 10 14 W/cm 2 Hence arbitrarily chosen ini- 
tial conditions in the full space will typically not lead to 
double ionization and numerical simulations of the whole 
process are rather unattainable. 

We have discussed the non-sequential double escape 
of the electrons considering trajectories within the sym- 
metry subspace. Motion in the symmetry subspace is 
unstable, that is deviations from the subspace will be 
amplified leading to single rather than double ionization. 
We can illustrate this with trajectories started slightly 
outside the symmetry plane (Fig. ^|). Fig ^|a shows ini- 
tial conditions on the saddle and symmetrically escaping 
electrons. For some deviation from symmetry, one elec- 
tron escapes, the other remains trapped to the nucleus 
(Fig. |b). 

It is possible, however, that the second electron returns 
to the nucleus but picks up enough energy to ionize itself 
(Fig. ^|c). In the adiabatic picture, if the energy of the 
remaining electron is higher than the saddle for a sin- 
gle electron atom, Vf — — 2y / 2\F\ 1 whether the electron 
stays trapped or escapes depends on the other constants 
of motion (besides the energy) in this integrable system. 
The second electron can thus escape during the same half 
field cycle as the first one but its final parallel momen- 
tum component need not be related to that of the first 
electron, as is shown in Fig. ^|c. 



100 



50 



-50 



(a) 



-100 



-50 50 100 150 



3 

d. o 

X -2 
-4 



(b) 



-10 



10 20 



20 
10 


-10 
-20 



(C) 




-40 -20 20 40 

z(a.u.) 

FIG. 9. Trajectories of electrons outside the symmetric 
subspace for E = —0.58 a.u. and F = 0.137 a.u.. Initial po- 
sitions are close to the saddle and marked by heavy dots; the 
electrons are distinguished by dotted and continuous tracks. 
Panel (a) shows a symmetric escape of the electrons. Panel 
(b) shows a case where outside the symmetry subspace one 
electron escapes and the other falls back to the ion. Panel (c) 
shows an example of sequential ionization of both electrons 
in opposite directions. 
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FIG. 10. Final parallel ion momenta distributions calcu- 
lated in the two-dimensional non-interacting electrons model 
for the initial time to = 0.33Ttj in the pulse duration, where 
T d — 20 x 2-k/lo, peak field amplitude F = 0.137 a.u. and 
initial energy E = —0.8 a.u. [panel (a)] and E = —0.1 a.u. 
[panel (b)]. The initial conditions for the electrons have been 
chosen to satisfy E = E± + E2 with the restrictions E\ < 
and E2 < 0. The saddle energy for a single electron is 
Vf = —1.05 a.u.. The results in each panel are based on 
about 1.5 • 10 6 trajectories. 

We are not able to calculate the contribution of sequen- 
tial events to the double ionization process but we can 
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estimate the distribution of the final ion momenta if the 
sequential process were dominant. To simulate sequen- 
tial escape one may return to a non-interacting electron 
model, essentially since the electrons cross the barrier at 
different times. In the model the initial conditions for 
each electron are chosen independently, constrained only 
by the requirement that the total initial energy is fixed 
E = Ei + E 2 and that E x < and E 2 < 0. Simulta- 
neous double ionization are not explicitly excluded, but 
events with delayed ionization are more probable, so that 
the distributions can still reflect the contributions from 
sequential ionization. 

In Fig. [l^ the distributions of parallel ion momenta 
for E = —0.8 a.u. and E = —0.1 a.u., calculated in the 
non-interacting electrons model for F = 0.137 a.u, are 
plotted. The figure should be compared with the figures 
from the previous section. The conclusion is straightfor- 
ward: the non-interacting electrons model is not able to 
reproduce correlations between the electrons observed in 
the experiments. Moreover, sequential escape can not be 
a dominant mechanism for double ionization in the range 
of the field intensities considered here. 



VII. CONCLUSIONS 

We have considered the process of double ionization 
of atoms in a strong, linearly polarized field for inten- 
sities below the saturation of single electron ionization. 
We have developed a model for non-sequential double 
ionization within classical mechanics. The process has 
been divided into two stages: in the first one a rescat- 
tering process leads to a highly excited complex of two 
electrons. In the second stage, an ionization of such com- 
pound state takes place. We have focused on the latter 
stage considering different pathways to double escape of 
the electrons. 

The excited complex can doubly ionize even when its 
energy is negative because the external field opens up 
saddles for electron escape. The pathway favored by 
the Coulomb interactions and the field is simultaneous 
symmetric escape of both electrons. Deviations form the 
symmetric configurations are amplified by the repulsion 
between the electrons which pulls one electron to infinity 
but the other one is pushed back to the nucleus. There- 
fore we propose that non-sequential double ionization is 
dominated by motions of electrons in or near the sym- 
metric subspace with the saddle. 

The requirement of the symmetric motion greatly sim- 
plifies the analysis which then can be carried out for the 
three- or even two-dimensional effective potential. The 
trajectory simulations within the symmetric configura- 
tions turns out to reproduce the experimentally observed 
ion momenta distributions very well. We have also con- 
sidered an alternative mechanism of the ionization, i.e. 
sequential escape of the electrons. By means of the non- 
interacting electrons model we show that the sequential 



ionization is not able to explain the experimentally ob- 
served electrons correlations. 

The modeling of the experimental distributions re- 
quires information on two parameters, the initial energy 
and the time of formation, which reflect a lack of knowl- 
edge on the compound state and the ramping of the field. 
The dependence of the momentum distributions on the 
parameters and comparison with the experimental results 
give insights into the dynamics of double ionization. 

The analysis in the present paper has been restricted 
to double ionization but its extension to multiple escape 
is straightforward p9| . In its current form the model 
is applicable for linearly polarized fields only. For other 
polarizations the number of rescattering events is greatly 
reduced. However, for some elliptically polarized field, 
if an electron is driven back to the core and a highly 
excited complex is formed, in the adiabatic approxima- 
tion the symmetric configuration of the electrons can be 
defined with respect to the temporary electric field axis. 
Then one can proceed with the analysis as for the linearly 
polarized case. 

Our whole discussion has a more than superficial sim- 
ilarity with Wanniers analysis of double ionization by 
electron impact II3-|l3]. The main differences are the 
presence of a field and its time dependence, which en- 
larges phase space and complicates the identification of 
the transition state. In the adiabatic approximation at 
fixed field strength we could identify this saddle in the 
C 2v subspace. The comparison with experiments is com- 
plicated furthermore by the necessity to average over 
initial energy and time of preparation of the compound 
state. Thus, signatures one might attribute to Wanniers 
analysis, such as threshold exponents (they follow imme- 
diately from the stability analysis of the saddle, but are 
not easy to verify), will be even more difficult to study. 
But we have no doubt that the observation of the corre- 
lated escape of the electrons [|J is a clear signature of the 
existence and dominance of the saddle and the pathways 
to double ionization which we discuss here. 
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